Revisiting the Excitation Energy Transfer in the Fenna-Matthews-Olson Complex 
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It is believed that the quantum coherence itsell cannot explain the very high excitation energy 
transfer (EET) efficiency in the Fenna-Matthews-Olson (FMO) complex. In this paper, we show 
that this is not the case if the inter-site couplings take complex values. By phenomenologically 
introducing phases into the inter-site couplings, we obtain the EET efficiency as high as 0.8972 
in contrast to 0.6781 with real inter-site couplings. Dependence of the excitation energy transfer 
efficiency on the initial states is elaborated. Effects of fluctuations in the site energies and inter-site 
couplings are also examined. 
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In 1962, John Olson isolated a water-soluble bacteri- 
ochlorophyll (BChl a) protein from green sulfur bacteria 
[l[. In 1975, Roger Fenna and Brian Matthews resolved 
the X-ray structure of this protein (Fenna-Matthews- 
Olson (FMO)) from Prosthecochloris aestuarii and found 
that the protein consists of three identical subunits re- 
lated by C3 symmetry, each containing seven BChl a 
pigments • In photosynthetic membranes of green sul- 
fur bacteria, this protein channels the excitations from 
the chlorosomes to the reaction center. Since it was the 
first photosynthetic antenna complex of which the X-ray 
structure became available, it triggered a wide variety 
of studies of spectroscopic and theoretical nature, and it 
therefore has become one of the most widely studied and 
well-characterized pigment-protein complexes. 

The excitation-energy transfer can be elucidated with 
2D echo-spectroscopy using three laser pulses which hit 
the sample within several femtoseconds time spacings. 
The experimentally measured 2D echo-spectra of the 
complex show wave-like energy transfer with oscillation 
periods roughly consistent with the eigenenergy spac- 
ings of the excited system |3[. This has brought a long- 
standing question again into the scientific focus that 
whether nontrivial quantum coherence effects exist in 
natural biological systems under physiological conditions 

Many studies have attempted to unravel the precise 
role of quantum coherence in the EET of light-harvesting 
complexes [(14151]. The interplay of coherent dynamics, 
which leads to a derealization of an initial excitation 
arriving at the FMO complex from the antenna, and 
the coupling to a vibronic environment with slow and 
fast fluctuations, has led to studies of environmentally 
assisted transport in the FMO. An interesting finding 
is that the environmental decoherence and noise play a 
crucial role in the excitation energy transfer in the FMO 
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In these studies, the FMO complex is treated using the 
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where |j) represents the state with only the j-th site be- 
ing excited and all other sites in their electronic ground 
state. Ej is the on-site energy of site j, and Jjj denotes 
the inter-site coupling between sites i and j. The inter- 
site couplings Jij given by Jy = J dr<&i(f)pj(r) represent 
the Coulomb couplings between the transition densities 
of the BChls [20(, where $i(r) is the electrostatic poten- 
tial of the transition density of BChl i and Pj(f) is the 
transition density of BChl j. It is easy to find that this 
calculation can only give real inter-site couplings Jy. In 
fact, the couplings, which has been used in previous stud- 
ies of 2D echo-spectra and facilitate comparisons between 
different approaches, are all real. 

We should emphasize that the inter-site couplings orig- 
inate from dipole-dipole couplings between different chro- 
mophores. They are in general complex from the aspect 
of quantum mechanics, regardless it is difficult to cal- 
culate. In this paper, we will shed light on the effects 
of complex inter-site couplings by addressing the follow- 
ing questions. First, we present a theoretical framework 
which elucidates how complex couplings could assist the 
excitation transfer. We optimize the relative phases in 
the couplings for maximal transfer efficiency and study 
the robustness of excitation transfer against the fluctua- 
tion of the on-site energy and inter-site couplings. Sec- 
ond, we present a study on the coherence in the initially 
excited states by showing the dependence of the transfer 
efficiency on the initial states. 

This paper is organized as follows. In Sec. II, we 
introduce the model to describe the FMO complex for 
the dynamics of the exciton transport. In Sec. Ill, 
we optimize the phases phenomenologically introduced 
in the couplings for maximal excitation energy transfer 
efficiency and explore how the efficiency depends on the 
initial state of the FMO. The fluctuations in the local 



2 



on-site energy and in the inter-site couplings affect the 
transfer efficiency, this effect is studied in Sec. IV. Fi- 
nally, we conclude in Sec. V. 

II. MODEL 

Since the FMO complex is arranged as a trimer with 
three different subunits interacting weakly with each 
other, we restrict our study to a single subunit which con- 
tains eight bacteriochlorophyll molecules (BCfils) [2ll. l22jj . 
The presence of the 8th BChl chromophore has just been 
suggested by the recent crystallographic data[21[, and 
the recent experimental data and theoretical studies in- 
dicated that the 8th BChl is the closest site to the base- 
plate and should be the point at which energy flows into 
the FMO complex [Hllll. However, the eighth BChl is 
loosely bound, it usually detaches from the others when 
the system is isolated from its environment to perform ex- 
periments. Thus we do not take the eighth into account 
and the EET will be explored by using the seven-site 



model[H|. 

We describe the FMO complex by a generic one- 
dimensional Frenkel exciton model, consisting of a reg- 
ular chain of 7 optically active sites, which are modeled 
as two-level systems with parallel transition dipoles. The 
corresponding Hamiltonian reads, 

H=Y i E j \m+ E (Ji^ui+jmiii), (2) 

j=l i>j = l 

where \ j) represents the state where only the j-th site is 
excited and all other sites are in their electronic ground 
state. Ej is the on-site energy of site j, and Jij denotes 
the excitonic coupling between sites i and j, Jij = J* . 
The on-site and the inter-site energies from Ref.[5] (in 
units of cm -1 ) will be adapted to study the quantum 
dynamics of the FMO. As afromentioned, a phase fcj 
is added to the inter-site coupling Jjy stronger than 15 
cmT x , 
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Here the zero energy has been shifted by 12230 cm -1 for 
all sites, corresponding to a wavelength of ~ 800nm. We 
note that in units of Ti = 1, we have 1 ps _1 =5.3 cm" 1 . 
Then by dividing |J,-j| and Ej by 5.3, all elements of 
the Hamiltonian are rescaled by units of ps _1 . We can 
find from the Hamiltonian H that in the Fenna-Matthew- 
Olson complex (FMO), there are two dominating exciton 
energy transfer (EET) pathways: 1 — > 2 — > 3 and 6 — > 
(5,7) — > 4 — > 3. Although the nearest neighbor terms 
dominate the site to site coupling, significant hopping 
matrix elements exist between more distant sites. This 
indicates that coherent transport itself may not explain 
why the excitation energy transfer is so efficient. 

We phenomenologically introduce a decoherence 
scheme to describe the excitation population transfer 
from site 3 to the reaction center (site 8). The master 
equation to describe the dynamics of the FMO complex 
is then given by, 

^ = -i[H,p]+CM, (4) 

where C 38 (p) = T [P 83 p(t)P 3 8 - \Pzp(t) - ±p(i)P 3 ] with 
Pss = |3)(8| = Pg 3 characterizes the excitation trapping 



at site 8 via site 3. 

We shall use the population P 8 at time T in the reac- 
tion center given by P 8 (T) = Tr(|8)(8|p(T)) to quantify 
the excitation transfer efficiency. Clearly, the Liouvillian 
-^38 (p) plays an essential role in the excitation transfer. 
With this term, we will show in the next section that the 
complex inter-site coupling can enhance the excitation 
energy transfer. 



III. SIMULATION RESULTS 

It has been shown that the Hamiltonian with real inter- 
site couplings can not give high excitation energy trans- 
fer efficiency. Indeed, the previous numerical simula- 
tions shown that T can be optimized to 87.14 with all 
4>ij = to obtain the highest excitation transfer efficiency 
Ps = 0.6781 in this case |24| . 

With complex inter-site couplings, the EET efficiency 
can reach almost 90% by optimizing the phases added to 
Jij. In units of 7r, these phases are (/>i2 = 1.2566, </>i6 = 
3.3510,023 = 1.8850,034 = 1.8850,045 = 1.88 50, 4 7 = 
0.8378,056 = 3.1416,0 6 7 = 0.3142. With these phases, 
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FIG. 1: Population on each site as a function of time. Site 
8 represents the reaction center. The phases are optimized 
for the EET efficiency. The exciton is assumed initially on 
the site 1. The upper panel shows the population without 
decoherence, while the bottom panel shows that with 72 = 2.4 
and 7j^2 = 0, see Eq.([5]). V — 44 is chosen throughout this 
paper. 



the population on each site as a function of time is plot- 
ted in Fig. [T] Two observations can be made from Fig. 
[TJ (1) The population shows oscillatory behavior in the 
whole energy transfer process, this is different from the 
situation where decoherence dominates the mechanism of 
EET [24|], (2) the excitation on site 1 and 2 dominates the 
population, while the population at site 3 is almost zero. 
The first observation can be understood as the quantum 
coherence lasts the whole EET, and the second observa- 
tion suggests that the site 1 and site 2 play important 
role in the excitation energy transfer. The excitation on 
site 3 is trapped and transferred to the reaction center al- 
most immediately as to obtain a high transfer efficiency, 
as Fig. [1] shows. 

Spatial and temporal relaxation of exciton shows that 
the site 1 and 6 were populated initially with large con- 
tribution [f|. Then it is interesting to study how the 
initial states affect the excitation transfer efficiency. In 
the following, we shall shed light on this issue. Two cases 
are considered. First, we calculate the excitation trans- 
fer efficiency with exciton initially in a superposition of 
site 1 and 6 and analyze the effect of transfer efficiency 
on the initial states. Second, we extend this study to 
initial states where sites 1 and 2 are initially excited. 
These analyses are performed by numerically calculate 
the transfer efficiency at time T = 5ps, selected results 
are shown in Fig. [5] and Fig. [3] 




FIG. 2: The dependence of the transfer efficiency on the ini- 
tial states. Top panel: the initial state site is a superposition 
of |1> and |6). Namely, \ip(t = 0)) = cos 0|1) +sin 9 exp(i<j>)\&) . 
Lower panel: the initial state is \tp(t = 0)) = cos0|l) + 
sin#exp(i0)|2). The phases in Jij are optimized for maximal 
transfer efficiency, i.e., they take the same values as in Fig[T] 
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FIG. 3: The transfer efficiency versus initial states: the case of 
classical superposition. The phases (couplings) are the same 
as in Fig. Q] The initial state is p|l)(l| + (l-p)|6)(6| for solid 
line, and it is p]l)(l| + (1 - p)|2)(2] for dotted line. 

In Fig. [2] we present the transfer efficiency as a 
function of and 0, which characterize the pure ini- 
tial states of the FMO complex through \ip(t — 0)} = 
cos0|l) + sin0exp(i0)|6) (upper panel) and \ip(t = 0)) = 
cos0|l)+sin0exp(i0)|2) (lower panel) . We find the FMO 
complex efficient to transfer excitation initially at site 1 
is not good at EET with site 6 occupied. Namely, coher- 
ent superposition of sites 1 and 6 decreases the exciton 
transfer efficiency. For exciton initially in a superposi- 
tion of 1 and 6, the transfer efficiency is more sensitive 
to the population ratio (characterized by 0) but not to 
the relative phase 4>. For exciton initially excited on site 
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71 = 2.4.73 = 2.0,76 
yields P s = 0.9374. 



1.8, the other 7, = and T = 44 



FIG. 4: EET efficiency as a function of 72. The other deco- 
herence rates jj are zero, the phases in the inter-site couplings 
are optimized for EET efficiency with 72 = 0. 



1 and 2, both the population ratio and the relative phase 
affect the energy transfer, the transfer efficiency reaches 
r7r and 4> — \k and it arrives at 



its maximum with 9 



its minimum with = ? and 



i-7r. This observation 
suggests that a properly superposition of site 1 and site 2 
can enhance the EET, although the increase of efficiency 
is small. 

Fig. [3] shows the dependence of the transfer efficiency 
on the mixing rate p, where the initial state is a classical 
mixing of site 1 and 6 (or 2). Obviously, the population 
mixing of site 1 and site 6 does not favors the transfer 
efficiency. While the classical mixing of site 1 and site 2 
increase the transfer efficiency a little bit (not clear from 
the figure). 

It is recently suggested that decoherence is an essential 
feature of EET in the FMO complex, where decoherence 
arises from interactions with the protein cage, the reac- 
tion center and the surrounding environment. Further 
examination shown that dephasing can increase the ef- 
ficiency of transport in the FMO, whereas other types 
of decoherence (for example dissipation) block the EET. 
These effects of decoherence can be understood as the 
fluctuation-induced-broadening of energy levels, which 
bridges the on-site energy gap and changes equivalently 
the coupling between the sites. 

We now examine if the decoherence can further in- 
crease the EET efficiency with the complex inter-site cou- 
plings. The decoherence can be described by an added 
Liouvillian term 



7j 



(5) 



to Eq. © with Pj = j = 1,2,3, 7. The de- 

coherence may come from site-environment couplings, 
where the environment models the thermal phonon bath 
and radiation fields. Optimizing the decoherence rates 
7,-,(j = 1,2,..., 7) with r = 44 for the EET efficiency, 
we find 72 = 2.4 and jj = 0,(j ^ 2) yields a transfer 
efficiency P 8 = 0.9344, see Fig. Q] Non-zero Jj,(j = 
1,3,4,5,6,7) can increase the EET efficiency, but the 
amplitude of the increase is not evident, for example, 



IV. EFFECT OF FLUCTUATIONS 

The site energies are among the most debated prop- 
erties of the FMO complex. These values are needed 
for exciton calculations of the linear spectra and simu- 
lations of dynamics. They depend on local interactions 
between the BChl a molecule and the protein envelope 
and include electrostatic interactions and ligation. Since 
the interactions are difficult to identify and even harder 
to quantify, the site energies are usually treated as in- 
dependent parameters obtained from a simultaneous fit 
to several optical spectra. There are many approaches 
to obtain the site energies by fitting the spectra. One of 
the main differences between those approaches is whether 
they restrict the interactions to BChl a molecules within 
a subunits or whether they include interactions in the 
whole trimer. Since the site energies and inter-site cou- 
plings may be different and fluctuate due to the couplings 
to its surroundings, it is interesting to ask how the fluctu- 
ations in the site energies and couplings affect the trans- 
fer efficiency. In the following, we will study this issue 
and show that the exciton transfer in the FMO complex 
remains essentially unaffected in the presence of random 
variations in site energies and inter-site couplings. This 
strongly suggests that the experimental results recorded 
for samples at low temperature would also be observable 
at higher temperatures. 

To be specific, we consider two types of fluctuations 
in the site energies and inter-site couplings. In the first 
one, it has zero mean, while the another type of fluc- 
tuations has nonzero positive mean. For the fluctuations 
with zero mean, the Hamiltonian in Eq. ([3]) takes the fol- 
lowing changes, Hjj — > Hjj(l + r\ ■ (rand(l) — 0.5)) and 
Hi& -> Hift^v +r 2 ■ (rand(l) -0.5)). For the fluctua- 
tions with nonzero positive mean, the Hamiltonian takes 
the following changes, Hjj —> Hjj(l + 77 ■ rand(l)) and 
H i= £j — > Hi^jie 1 ^ +r 2 -rand(l)), where rand(l) denotes 
a random number between and 1. So a 100% static 
disorder may appear in the on-site energies and inter-site 
couplings. 

With these arrangements, we numerically calculated 
the transfer efficiency and present the results in Fig. [5] 
Each value of transfer efficiency is a result averaged over 
20 fluctuations. From the numerical simulations, we no- 
tice that the results for fluctuations with non-zero mean 
are almost the same as those for zero-mean fluctuations, 
we hence present here only the results for zero-mean fluc- 
tuations. Two observations are obvious. (1) With zero- 
mean fluctuations in the site energies and couplings in- 
crease, the transfer efficiency fluctuates greatly. (2) The 
efficiency decreases with both r-i and 77. Moreover, the 
EET efficiency is sensitive to the fluctuations in the site 
energy more than that in the inter-site couplings. 

These feature can be interpreted as follows. Since the 
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FIG. 5: Effect of fluctuations in the Hamiltonian on the trans- 
fer efficiency. The plot is for the fluctuations with zero mean. 
The resulting Ps is an averaged result over 20 independent 
runs. The fluctuations with non-zero mean has a similar ef- 
fect on the EET, which is not shown here. 



network with Hamiltonian Eq.© is optimal for EET 
efficiency, any changes in the site-energy and inter-site 
couplings would diminish the excitation energy transfer, 
leading to the decrease in the EET efficiency. This fea- 
ture is different from the case that the decoherence dom- 
inates the mechanism of EET. As Ref. [HJ shown, the 
efficiency increases with r\ but decreases with ti- The 
physics behind this difference is as follows. For the case 
that the decoherence dominates the mechanism of EET, 
the energy gap between neighboring sites blocks the en- 
ergy transfer, whereas the inter-site couplings (represent- 
ing the overlap of different sites) together with the de- 
coherence favor the transport. Whereas for the present 
case the phases in the inter-site couplings need to match 
the energy spacing between sites. As a result, any mis- 
matches would results in decrease in the EET efficiency. 

It was suggested that the fluctuations in the on- 
site energies are Gaussian. To examine the effect of 
Gaussian fluctuations, we introduce a Gaussian function, 

_(x- M ) 2 

y(x\a, fx) = ^J^ e 2 "' 2 i where fi is the mean, while a 
denotes the standard deviation. We assume that the fluc- 
tuations enter only into the on-site energies, namely, Hjj 
is replaced by Hjj(l + y(x\a, //)), j — 1,2, ...,7. With 
these notations, we plot the effect of fluctuations on the 
transfer efficiency in Fig. [Scupper panel). We find that 
the dependence of the efficiency on the variation is subtle: 
When the mean is zero, the efficiency reach its maximum 
at a = 0.5, while when the mean is large, small variation 
favors the EET efficiency. Moreover, we find that the ef- 
ficiency decreases as the mean increases. This again can 
be understood as mismatch between the site energy and 
the inter-site couplings caused by the fluctuations. When 
these Gaussian fluctuations occur in both the on-site 
energies and the couplings, we find from Figinjbottom 
panel) that the transfer efficiency decreases as the mean 
and variation increases. 

Absorption spectroscopy and 2D echo-spectroscopy are 
used to study time-resolved processes such as energy 
transfer or vibrational decay as well as to measure in- 




FIG. 6: Effect of fluctuations in the Hamiltonian on the trans- 
fer efficiency. The fluctuations are Gaussian with mean \i and 
the standard deviation a. Upper: The fluctuations happen 
only in the on-site energies; Bottom: The fluctuations in both 
the on-site energies and couplings. The resulting fidelity is an 
averaged result over 15 independent runs. 



termolecular couplings strengths. They provide a tool 
that gives direct insight into the energy states of a quan- 
tum system. In turn, the eigenenergies are essential 
to calculate the absorption spectroscopy and 2D echo- 
spectroscopy. It is easy to find that the eigenenergies 
and energy gaps for the FMO with complex inter-site 
couplings, given by (96.2851, 62.2958, 61.4307, 48.2046, 
22.6379, 18.9438, -4.1374), are almost the same as in the 
system with real inter-site couplings given by (96.8523, 
62.6424, 57.9484, 50.6357, 22.8218, 19.2387, -4.4789). 
The maximal difference is at the 7th eigenenergy, it is 
about 8.25% to its eigenenergy. This suggests that the 
absorption spectroscopy and 2D echo-spectroscopy may 
not change much due to the introduced phases in the 
inter-site couplings. 



V. CONCLUSION 

Excitation energy transfer (EET) has been an inter- 
esting subject for decades not only for its phenomenal 
efficiency but also for its fundamental role in Nature. 
Recent studies show that the quantum coherence itself 
cannot explain the very high excitation energy transfer 
efficiency in the Fenna-Matthews-Olson Complex. In this 
paper, we show that this is not the case when the inter- 
site couplings are complex. Based on the Frenkel exciton 
Hamiltonian and the phenomenologically introduced de- 
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coherence, we have optimized the phases in the inter-site 
couplings for maximal energy transfer efficiency, which 
can reach about 89.72% at time 5ps without decoherence 
and 93.44% with only dephasing at site 2. By considering 
different mixing of exciton on site 1 and site 6 (site 2) as 
the initial states, we have examined the effect of initial 
states on the energy transfer efficiency. The results sug- 
gest that any mixing of site 1 and site 6 decreases the 
energy transfer efficiency, whereas a coherent superposi- 
tion (or classically mixing) of site 1 and site 2 does not 
change much the EET efficiency. The fluctuations in the 
site energies and inter-site couplings diminish the EET 
due to the mismatch caused by the fluctuations. 

Finishing this work, we have noticed a closely related 



preprint . While they focus on bi-pathway EET in the 
FMO (i.e., neglecting the inter-site couplings weaker than 
15 cm -1 ) and a general representation for complex net- 
work, we are mainly concerned about the phases added 
to the inter-site couplings, initial excitations and fluctua- 
tions in the site energies and inter-site couplings. In this 
respect, both works are complementary to each other. 
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